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1. Introduction 

The most important discovery in particle physics in the past decade has only deepened the 
mystery of "flavor" of quarks and leptons. The fact that the mixing angles in the leptonic 
sector are large ]l], |2| stands in sharp contrast with the observed small mixing angles in the 
quark sector. In particular, mixing between the second and third generation neutrinos appears 
to be maximal. Of course, this large mixing could occur from diagonalizing the neutrino mass 
matrix, the charged lepton mass matrix, or both. At present, the source of this large mixing 
is a mystery. 

In view of this, it is tempting to explore other interactions which change lepton flavor 
between the second and third generations. Several years ago, two of us (TH, MS), along with 
Black and He (BHHS) ||, performed a comprehensive analysis of constraints on these inter- 
actions based on low energy meson physics. BHHS chose an effective field theory approach, 
in which all dimension-6 operators of the form 

GurT)(«fiY), (i.i) 
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were studied, where F contains possible Dirac 7-matrices. With six flavors of quarks, there 
were 12 possible combinations of q a and q b (assuming Hermiticity) , six diagonal and six off- 
diagonal, and four choices S, P, V, A of the gamma matrices were considered. All of these 
operators were considered, and most were bounded by considering r, K, B and t decays. 
In particular, BHHS considered operators of the form 

AC = A45 = J2^-fc Ii r) (r + H.c. , (1.2) 

where Tj £ (1, 7 5 , 7 CT , 7 CT 7 5 ) denotes relevant Dirac matrices, specifying scalar, pseudoscalar, 
vector and axial vector couplings, respectively. They did not consider tensor operators since 
the hadronic matrix elements were not known and the bounds were expected to be weak in 
any event. They chose a value of 

C{ p = 4tt0(1) (default), (1.3) 

which corresponds to an underlying theory with a strong gauge coupling of as = 0(1). 
Arguments can be made for multiplying or dividing this by 4tt, for naive dimensional analysis 
or for weakly coupled theories, respectively. A discussion is found in BHHS; we simply choose 
the above definition of A and other choices can be made by simple rescaling. 

Besides the four fermion operators in Eq. ( |1.2j ), there may be other induced operators 
involving the SM gauge bosons, such as the electroweak transition operator 

AC = ^Jut^tF^, (1.4) 

where v is the vacuum expectation value of the Standard Model Higgs field and Fn V is the 
electroweak field tensor. However, when these operators are compared to the underlying 



new strong dynamics of the four fermion interaction in Eq. (1.2), it is found that they are 
suppressed by 0(Mw/A), where M\y is the mass of the electroweak gauge boson. For new 
physics scales of order 1 TeV or greater, this is at least an order of magnitude suppression. 
Thus, we ignore these operators. 

BHHS found that operators involving the three lightest quarks were strongly bounded, 
with bounds ranging from 3 to 13 TeV on the related value of A. These bounds can be found 
in Appendix [A|. Not surprisingly, operators involving the top quark were either unbounded or 
very weakly bounded, with only the tu operator for vector and axial vector couplings being 
bounded by A < 650 GeV (the bound arises through a loop in B — > \it decay). Operators 
involving the 6-quark and a light quark also have bounds on A which were generally in 
the several TeV range. However, there were some surprises. The scalar and pseudoscalar 
operators involving cu and cc were completely unbounded, and the bb operator was essentially 
unbounded for all S, P, V, A operators. And, as noted above, none of the tensor operators 
were considered at all, for all quark combinations. 

In this note, we point out that the operators in Eq. ( |1 . 1|) (without involving top quarks) 
will contribute to fj, — r production at hadron colliders. Given that many of the possible 
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operators, as noted above, are completely unbounded or weakly bounded from the current 
low energy data, study of pp — > fir at the LHC or pp — > [it at the Tevatron will probe 
unexplored territory. 

There have been some previous discussions of /i — r production at hadron colliders. Han 
and Marfatia [|j looked at the lepton- violating decay h — > /ir at hadron colliders, and a very 
detailed analysis of signals and backgrounds was carried out by Assamagan et al. ||] after- 
wards. Other work looking at Higgs decays focused on mirror fermions ||, supersymmetric 
models pj, seesaw neutrino models (8|, and Randall-Sundrum models H. In addition to 
Higgs decays, others have considered lepton-flavor violation in the decays of supersymmetric 



particles and in horizontal gauge boson models [11]. These analyses, however, were done 
in the context of very specific models (often relying on the assumption that the /i and r are 
emitted in the decay of a single particle). Here, we will use a much more general effective 
field theory approach. 

This paper is organized as follows. In the next section, we discuss the cross sections 
for /it production via the various operators. A detailed analysis of the signal identification 
and background subtraction is in Section 3, and Section 4 contains some discussions and our 
conclusions. Appendix [A] reiterates the bounds from BHHS for comparison, and Appendix [B| 
outlines the calculation of partial-wave unitarity bounds. 

2. /ir Production at Hadron Colliders 

Due to the absence of appreciable [it production in the SM, their production can be estimated 
via the effective operators in Eq. (|1,1|). On dimensional grounds, the cross section for qiqj — > 
\xt grows with center of mass energy, i.e., 

viQiQj -»• l*r) cx -p, (2.1) 

where y/s is the center of mass energy for the partonic system. This growth of cross section 
with energy will eventually violate unitarity bounds. Expanding the scattering amplitudes in 
partial waves, we find the unitarity bounds to be (see Appendix |b|) 

I 2A 2 for scalar, pseudoscalar, and tensor; 

S < < n (2.2) 

I 3A vector and axial vector case. 

The total cross sections for fir production at the hadronic level after convoluting with 
the parton distribution functions (pdfs) are 

^Scalar = J d,T (q <8> ?) (t) ^1 - — J T (2.3) 

Sector = dT{q ® q){T) { l ~T)\ l + ^) T ^ 

^Tensor = y f"" dr(q ® q) (r) (l - ^ (l + 2^ r, (2.5) 
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where r = s/S, to = rn^/S, m T is the tau mass, and \/S is the center of mass energy in the 
lab frame. The pseudoscalar cross section is of the same form as the scalar cross section, and 
the axial vector cross section is of the same form as the vector cross section. Our perturbative 
calculation will become invalid at the unitarity bound, hence there is a maximum on the r 
integration. It is given by r max = 2A 2 / S for the scalar, pseudoscalar, and tensor cases, and 
Tmax = 3A 2 /S for the vector and axial- vector cases. Also, q(x) is the quark distribution 
function with flavor sum suppressed, and ® denotes the convolution defined as 

(91 ®92)(y) = / dxi / dx 2 gi(xi)g 2 (x 2 )S(xix 2 - y). (2.6) 



J 



The CTEQ6L parton distribution function set is used for all of the results [12|. 

Results for the cross sections for the scalar, pseudoscalar, vector, axial vector, and tensor 
structures at the Tevatron, LHC at 10 TeV and 14 TeV are given in Table |l|. The cross section 
for the pseudoscalar (axial vector) current is the same as for the scalar (vector) current. For 
all cases, A is set equal to 2 TeV and the unitarity bounds are taken into consideration. At 
this rather high scale, the production rates are dominated by the valence quark contributions. 
The cross sections at the LHC are larger than those at the Tevatron by roughly an order of 
magnitude, reaching about 100 pb. 

For some cases the bounds from BHHS are greater than 2 TeV, hence the cross section 
needs to be scaled to determine a realistic cross section at hadron colliders. The partonic 
cross sections scale at A -4 , but at the hadronic level a complication arises since the unitarity 
bounds introduce a dependence on the new physics scale in the integration over pdfs. If 
the unitarity bounds are ignored (T max = 1), one finds that with A = 2 TeV neglecting the 
unitarity bounds has at most a 10% effect on the cross sections at the LHC for both 10 TeV 
and 14 TeV and no effect at the Tevatron since the unitarity bounds are greater than the lab 
frame energy. Hence, if A is increased from 2 TeV, at the LHC it is a good approximation to 
assume the cross section scales as A -4 and at the Tevatron the cross section scales exactly as 
A -4 . For example, the lower bound on A for the vector uu coupling from BHHS is 12 TeV, so 
the maximum cross section at the 14 TeV LHC from this operator would be approximately 
160 x (2/12) 4 pb = 120 fb. On the other hand, there is no bound whatsoever for the vector 
uc coupling, and thus a cross section limit of 110 pb would yield a new limit of 2 TeV on the 
scale of this operator. This would constitute an improvement of many orders of magnitude. 



3. Signal Identification and Backgrounds 

Upon production at hadron colliders, r's will promptly decay and are detected via their decay 
products. About 35% of the time the r decays to two neutrinos and an electron or muon, the 
other 65% of the time the r decays to a few hadrons plus a neutrino. We will consider the r 
decay to an electron as well as hadronic decays in this work. The decay to a muon will result 
in a final state that has a large Drell-Yan background. We will study the signal reach 

at the Tevatron and at the 14 TeV LHC. 
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Table 1: Cross sections for all the scalar, pseudoscalar, vector, axial vector, and tensor structures at 
the Tevatron at 2 TeV, the LHC at 10 TeV, and the LHC at 14 TeV. The pseudoscalar (axial vector) 
cross section is the same as the scalar (vector) cross section. All cross sections were evaluated with 
the new physics scale A = 2 TeV and the unitarity bounds are taken into consideration. 





Tevatron 2 TeV (pp) 


LHC 10 TeV (pp) 


LHC 14 TeV (pp) 


a (pb) 


1,75 


7^,7^75 




1,75 


7^,7/^75 




1,75 


7/^,7^75 




uu 


8.4 


11 


22 


63 


85 


170 


120 


160 


310 


dd 


2.5 


3.3 


6.7 


38 


51 


100 


72 


98 


190 


ss 


0.18 


0.24 


0.49 


5.5 


7.4 


15 


11 


15 


30 


ds 


1.3 


1.7 


3.4 


34 


45 


91 


66 


89 


180 


db 


0.50 


0.67 


1.3 


17 


22 


45 


34 


46 


90 


sb 


0.13 


0.17 


0.34 


5.0 


6.7 


13 


11 


14 


28 


uc 


1.5 


2.0 


3.9 


41 


55 


110 


80 


110 


210 


cc 


0.070 


0.094 


0.19 


2.6 


3.5 


7.0 


5.5 


7.3 


15 


bb 


0.021 


0.028 


0.056 


1.1 


1.5 


2.9 


2.4 


3.2 


6.4 



3.1 r Decay to Electrons 
3.1.1 Signal Reconstruction 

The r decays to an electron plus two neutrinos about 18% of the time. We thus search for a 
final state of an electron and muon 

e + £t. (3.1) 

The electromagnetic calorimeter resolution is simulated by smearing the electron energies 
according to a Gaussian distribution with a resolution parameterized by 

a(E) a , . 

where the constants are a = 10% and b = 0% at the Tevatron |13|], a = 5% and b = 0.55% at 
the LHC [0, and ® indicates addition in quadrature. For simplicity, we have used the same 
form of smearing for the muons. 

The decay of the r leaves us with some missing energy and we need to consider how to 
effectively reconstruct the r momentum. For our process all the missing transverse momentum 
is coming from the r, hence 

„t „e i „miss /o o\ 

Pt = Pt + Pt ■ 

At hadron colliders, we have no information on the longitudinal component of the missing 
momentum on an event-by-event basis. However, the r will be highly boosted and its decay 
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products will be collimated. Hence, the missing momentum should be aligned with the 
electron momentum and the ratio p e z /p™ ms should be the same as the ratio of the magnitudes 
of the transverse momenta, pfp/p^ lss . Therefore, the longitudinal component of the r can be 
reconstructed as |Q] 

( miss \ 
1 + %r) • (3 - 4) 

Once the three- momentum is reconstructed, we can solve for the r energy, E% = + m^. 
Figure |l] illustrates the effectiveness of this method at the Tevatron. Figure 1(a) (Figure 1(b) ) 



shows the transverse momentum (longitudinal momentum) distribution for the theoretically 
generated (solid) and kinematically reconstructed (dashed) r momenta. As can be seen, the 
r momentum is reconstructed effectively. 

We first apply some basic cuts on the transverse momentum and the pseudo rapidity 
to simulate the detector acceptance and triggering, as well as to isolate the signal from the 
background, 

p% > 20 GeV, \r)»\ < 2.5, 

p e T > 20 GeV, \rf \ < 2.5. (3.5) 

Since the signal does not contain any jets, we also require a jet veto such that there are no 
jets with pt > 50 GeV and |r/| < 2.5. 

There are several distinctive kinematic features of our signal. The decay products of the 
r will be highly collimated, and the electron transverse momentum will be traveling in the 
same direction as the missing transverse momentum. Also, in the transverse plane the muon 
and tau should be back to back. Since the electron will mostly be in the direction of the r, 
it will also be nearly back to back with the muon. Finally, the r and \x have equal transverse 
momenta; hence, the decay products of the r have less transverse momentum than the [i. We 
can measure this discrepancy using the momentum imbalance 

A PT =p»-p e T . (3.6) 

For the signal, this observable should be positive. Based on the kinematics of our signal, we 
apply the further cuts f§ 

H(Pt,Pt) > 2-75 rad, 5<f){pf ss ,p e T ) < 0.6 rad, (3.7) 
Ap T > 0. 



3.1.2 Backgrounds and their Suppression 

The leading backgrounds are W + W~ pair production, Z^/j* — > t + t~ , and tt pair production 
H . The total rates for these backgrounds at the Tevatron and the LHC are given in Table § 
with consecutive cuts. We consider both of the final states with (i + and \x~ . 
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Figure 1: Distributions of the theoretically generated (solid line) and kinematically reconstructed 
(dashed line) r momentum at the Tevatron at 2 TeV with a uc initial state, scalar coupling, and new 
physics scale of 1 TeV. Fig. (a) is the r transverse momentum distribution, and Fig. (b) is the r 
longitidunal momentum distribution. 



Table 2: Leading backgrounds to the r's electronic decay before and after consecutive kinematic and 
invariant mass cuts for (a) the Tevatron at 2 TeV and (b) the LHC at 14 TeV. 



Backgrounds (pb) 


No Cuts 


Cuts Eq. 03.50 


+ Eq. @ 


+ Eq. © 


(a) Tevatron 2 TeV 


W + W~ -> n ± v il j^v T 


0.032 


0.0046 


0.0012 


2.6 x 10" 4 


W + W~ -»■ fl^U^Ve 


0.18 


0.13 


0.0060 


9.8 x 10" 4 


Z°/7* -»■ r+r" -»■ M ± ^ rT 


610 


0.21 


0.091 


1.4 x 10" 4 


— )• ^VybT^Vrb 


0.020 


6.5 x 10" 4 


7.4 x 10" 5 


4.4 x 10~ 5 


ti — > n Vftbe^vJ) 


0.11 


0.0099 


7.3 x 10" 4 


2.7 x 10" 4 


(b) LI 


iC 14 TeV 




0.34 


0.030 


0.0088 


0.0031 




1.9 


0.99 


0.051 


0.014 




2300 


1.1 


0.49 


0.0014 


— )• / u ± t' M 6r =F ^ T 6 


1.9 


0.070 


0.010 


0.0077 


tt — > ^v^be^vj) 


11 


1.5 


0.10 


0.050 



The partonic cross section of our signal increases with energy while the cross sections 
of the backgrounds will decrease with energy. Hence, the invariant mass distribution of our 
signal does not fall off as quickly as the backgrounds. 

Figure [2(a) shows the invariant mass distributions of backgrounds and our signal at the 
Tevatron with initial states cc and uc with various couplings and a new physics scale of 
1 TeV after applying the cuts in Eqs. ([^) and (|3.7|) . The cross section for the pseudoscalar 
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M (GeV) M (GeV) 

(XT v ' fIT v ' 

(a) (b) 

Figure 2: The invariant mass distributions of the reconstructed r — /i system at the Tevatron at 2 
TeV. Fig. (a) shows the distributions of the leading backgrounds (dotted and dot-dot-dash) and of 
our signal for the uc and cc initial states with coupling of various Lorentz structures and a new physics 
scale of 1 TeV. Fig. (b) shows the distributions of the leading backgrounds (dotted and dashed) and 
of our signal (solid) for the uc initial state with scalar coupling and various new physics scales. The 
cuts in Eqs. (3.5) and (3.7) have been applied. 



(axial-vector) couplings are the same as those for the scalar (vector) couplings. The decline 
in the signal rates is due to a suppression of the pdfs at large x. Although the signal rates 
steeply decline with invariant mass the background falls off faster. The uc signal is still clearly 
above background due to a valence quark in the initial state, but the cc signal distribution is 
much closer to the background distribution due to the steep fall with invariant mass and a 



lack of an initial state valence quark. Figure 2(b) shows the invariant mass distributions of 
backgrounds and our signal at the Tevatron with initial state uc and scalar coupling for various 
new physics scales. The 3 TeV new physics scale invariant mass distribution is approaching 
the background distribution. A higher cutoff on the invariant mass will be needed to separate 
the weak signal from the backgrounds. Based on Fig. H, we propose a selection cut on 



> 250 GeV. 



(3.8) 



Table shows the effects of the invariant mass cut on the backgrounds in the last column. 



Similar analyses can be carried out for the LHC. Figure |3(a)| shows the invariant mass 
distribution for our signal with the uc and cc initial states and various Lorentz structures, as 



well as the backgrounds after the cuts in Eqs. (3.5) and (|3.7|). The new physics scale was set to 



1 TeV and the unitarity bound is imposed. Figure 3(b) shows the invariant mass distribution 
of the uc initial state with various new physics scales. The cutoff on the invariant mass 
corresponds to the unitarity bound, the scale at which the perturbative calculation becomes 
untrustworthy. In the lack of the knowledge for the new physics to show up at the scale A, 
we simply impose a sharp cutoff at the unitarity bound. As compared with the Tevatron, 
the LHC signal rates fall off much less quickly with invariant mass since the Tevatron's lower 
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M (GeV) M (TeV) 

(XT (IT 

(a) (b) 

Figure 3: The invariant mass distributions of the reconstructed r — fx system at the LHC at 14 TeV. 
Fig. (a) shows the distributions of the leading backgrounds (dotted and dot-dot-dash) and of our 
signal for the uc and cc initial states with coupling of various Lorentz structures and a new physics 
scale of 1 TeV. Fig. (b) shows the distributions of the leading backgrounds (dotted and dashed) and 
of our signal (solid) for the uc initial state with scalar coupling and various new physics scales. The 
cut offs in the distributions at high invariant mass are due to the unitarity bounds. The cuts in Eqs. 
(3.5) and (3.7) have been applied. 



energy leads to a suppression from the pdfs at large x. As can be seen, as the new physics 
scale increases the cross section decreases and the background becomes more problematic at 
lower invariant mass. Also, as the new physics scale increases the unitarity bound becomes 
less strict. Hence, although the backgrounds at the LHC are considerably larger than at the 
Tevatron, for large new physics scales the LHC has an enhancement in the signal cross section 
from the large invariant mass region. 

3.2 r Decay to Hadrons 

Although with significantly larger backgrounds, the signal from r hadronic decays can be 
very distinctive as well. We limit the hadronic r decays to 1-prong decays to pions, i.e., 
r ± — > tt^Vt, — > Tr^~ir®i/ T , and — > tx^2t^v t . The r's have 1-prong decays to these final 
states about 50% of the time. We thus search for a final state of a r jet and a muon 

Jt + V- (3-9) 
To simulate detector resolution effects, the energy is smeared according to Eq. (ET^) with 



a = 80% and b = 0% for the jet at the Tevatron U and a = 100% and b = 5% at the LHC 
As in the electronic decay, the r is highly boosted and its decay products are collimated. 
Hence, all the missing energy in the event should be aligned with the r. The signal is then 
reconstructed as described in Eqs. ( |3.3|) and ( |3.4|) with the electron momentum replaced by 
the momentum of the r-jet. 
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The hadronic decay of the r also has the backgrounds W + W~ pair production, Z° /j* — > 
t + t~, and tt pair production plus an additional background of VK+jet, where the jet is 
misidentified as a r-jet. At the Tevatron, we assume a r-jet tagging efficiency of 67% and 
that a light jet is mistagged as a r-jet 1.1% of the time [|l^] and at the LHC we assume a r-jet 



tagging efficiency of 40% and a light jet misidentification rate of 1% [14|. Even with a low 
rate of misidentification, the VF+jet background is large. To suppress this background, we 
note that for hadronic decays most of the r transverse momentum will be carried by the jet. 
Hence the r-jet should be traveling in the same direction as the reconstructed r momentum. 



Motivated by this observation, we apply the same cuts as Eqs. (|3.5|) , (|3.7|) , and (|3.q) with the 
electron momentum replaced by the r-jet momentum and the additional cuts 

r—jet 

> 0.6 AR( P T~ jet ,Pt) < 0.2 rad. (3.10) 

Pt 

3.3 Sensitivity Reach at the Tevatron 

One can determine the sensitivity of the Tevatron to the new physics scale with 8 fb _1 of 
data. Table |3| shows the sensitivity of the Tevatron for (a) electronic and (b) hadronic r 
decays. The tables list the maximum new physics scale sensitivity at 2a and 5a level at the 
Tevatron. The reaches for scalar (vector) and pseudoscalar (axial-vector) are the same at the 
Tevatron, although the previous bounds from BHHS for the scalar (vector) and pseudoscalar 
(axial-vector) couplings may not be the same. The bounds from BHHS can be found in 
Appendix [A]. If only one of the bounds for scalar (vector) or pseudoscalar (axial-vector) 
coupling from BHHS is greater than the Tevatron reach one star is placed next to the new 
physics scale, if both bounds are greater than the Tevatron reach two stars are placed next 
to the new physics scale. Due to the larger backgrounds from W+jet, the Tevatron is much 
less sensitive to the r hadronic decays than the r electronic decays. 

There were no bounds from BHHS for the tensor couplings, so the Tevatron will be 
able to exlude some of the parameter space. Since the tensor cross sections are generally at 
least twice as large as the scalar cross sections, the Tevatron is more sensitive to the tensor 
couplings than it is to scalar couplings. Also, in general, the Tevatron is more sensitive to 
processes with initial state valence quarks than those without initial state quarks. With 8 
fb _1 of data most of the bounds can be increased, some quite stringently. 

Somewhat similar leptonic final states have been searched for in a model-independent 
way at the Tevatron [16|, although these included substantial missing energy and possible 



jets. We encourage the Tevatron experimenters to carry out the analyses as suggested in this 
article. 

3.4 Sensitivity Reach at the LHC 

The LHC is also sensitive to flavor changing operators. For the signal and background anal- 



ysis, we used the same kinematical cuts as we used at the Tevatron, see Eqs. ( |3.5[) , (|3.7|) , 



and (|3.8| ). Table |I] shows the sensitivity of the LHC to all possible initial states and the 
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Table 3: Maximum new physics scales the Tevatron is sensitive to with 8 fb _1 of data at the 2a 
and 5tr levels. The sensitivities are presented for both (a) electronic and (b) hadronic r decays with 
various initial states. One star indicates that the Tevatron reach is less than only one of the scalar 
(vector) or pseudoscalar (axial- vector) bounds from BHHS, and two stars indicates that the Tevatron 
reach is less than both bounds from BHHS. BHHS does not contain bounds on the tensor coupling. 



(a) r — > e 


Anp (TeV) 


2a sensitivity 


5cr discovery 


Coupling 


1,75 


7^,7^75 




1,75 


7/i,7/i75 




uu 


20 


21 


24 


14 


15 


17 


dd 


17 


18 


21 


12 


13 


15 


ss 


9.9 


10 


12 


7.2* 




8.7 


ds 


15 


16 


18 


10 


11* 


13 


db 


13 


14 


16 


9.8 


10* 


11 


sb 


9.5 


10 


11 


6.9 


7.3 


8.3 


uc 


17 


18 


20 


12 


13 


14 


cc 


7.9 


8.3 


9.5 


5.7 


6.0 


6.9 


bb 


6.4 


6.8 


7.7 


4.6 


4.9 


5.6 


(b) t — y 


Anp (TeV) 


2a sensitivity 


5a discovery 


Coupling 


1,75 


7^,7^75 




1,75 


7^,7^75 


CV 


uu 


8.6** 


9.2** 


10 


6.5** 


6.9** 


8.1 


dd 


5.7** 


6.1** 


7.1 


4.3** 


4.6** 


5.4 


ss 


1.8* 


1.9** 


2.3 


1.4** 


1.4** 


1.7 


ds 


3.7 


4.0* 


4.6 


2.8* 


3.0** 


3.5 


db 


2.7* 


2.9* 


3.4 


2.0** 


2.2 


2.5 


sb 


1.5** 


1.6** 


1.9 




1.2** 


1.4 


uc 


3.9 


4.1 


4.8 


2.9 


3.1 


3.6 


cc 


1.1 


1.2 


1.4 


0.89 


0.95** 


1.1 


bb 


0.91 


0.97 


1.1 


0.68 


0.73 


0.86 



couplings under consideration with 100 fb of data. The table contains the maximum new 
physics scales the LHC is sensitive to at the 2a and 5a levels. As with the Tevatron, the 
LHC reach for scalar (vector) couplings is the same as that for pseudoscalar (axial-vector) 
couplings, although the bounds from BHHS may be different. If only one of the bounds for 
scalar (vector) or pseudoscalar (axial-vector) coupling from BHHS is greater than the LHC 
reach one star is placed next to the new physics scale, if both bounds are greater than the 
LHC reach two stars are placed next to the new physics scale. Despite the larger backrounds 
for the hadronic r decays, at the LHC the reaches for the hadronic and electronic r decays are 
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Table 4: Maximum new physics scales the LHC is sensitive to at 14 TeV with 100 fb _1 of data at the 
2<T and 5cr levels. The sensitivities are presented for both (a) electronic and (b) hadronic r decays with 
various initial states. One star indicates that the LHC reach is less than only one of the scalar (vector) 
or pseudoscalar (axial-vector) bounds from BHHS, and two stars indicates that the LHC reach is less 
than both bounds from BHHS. BHHS does not contain bounds on the tensor coupling. 



(a) r — > e 


Anp (TeV) 


2a sensitivity 


5a discovery 


Coupling 


1,75 


7/^75 




1,75 


7^,7^75 




uu 


18 


19 


21 


14 


15 


17 


dd 


16 


17 


19 


12 


13 


15 


ss 


9.0* 


9.6* 


11 


7.1* 


7.6** 


8.6 


ds 


13 


14 


16 


10 


11* 


13 


db 


12 


13 


14 


9.7 


10 


11 


sb 


8.7 


9.2 


10 


6.8 


7.3 


8.2 


uc 


15 


16 


18 


12 


13 


14 


cc 


7.2 


7.6 


8.6 


5.7 


6.0 


6.8 


bb 


5.8 


6.2 


7.0 


4.6 


4.9 


5.5 


(b) r -> h ± 


Anp (TeV) 


2a sensitivity 


5a discovery 


uu 


15 


16 


18 


12 


13 


14 


dd 


13 


14 


16 


10* 


11* 


13 


ss 


7.9* 


8.4** 


9.7 


6.2* 


6.7** 


7.7 


ds 


11 


12* 


14 


9.3 


9.9* 


11 


db 


10 


11 


13 


8.4* 


8.9 


10 


sb 


7.6 


8.1 


9.3 


6.0 


6.4 


7.4 


uc 


13 


14 


16 


10 


11 


12 


cc 


6.3 


6.7 


7.8 


5.0 


5.3 


6.2 


bb 


5.1 


5.5 


6.3 


4.1 


4.3 


5.0 



much more similar than at the Tevatron since the LHC cross section receives an enhancement 
from the large invariant mass region. For electronic (hadronic) r decays the LHC with 100 
fb _1 of data is less (more) sensitive than the Tevatron with 8 fb _1 of data. 

Figure Q shows the integrated luminosities needed for 2a and 5a observation at the LHC 
with various initial states and r decay to electrons as a function of the new physics scale. 
For some initial states and Lorentz structures BHHS had a bound on the new physics scale 
larger than 1 TeV. In those cases the distribution does not begin until the BHHS bound on 
the new physics scale. The sensitivity for the pseudoscalar (axial-vector) is the same as the 
scalar (vector) state, although the bounds from BHHS are different. Note that extraordinary 
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Figure 4: The luminosity at the 14 TeV LHC needed for 2cr and 5cr observation as a function of the 
new physics scales with couplings of various Lorentz structures and electronic r decay. The sensitivity 
for the uc initial state is shown in (a), for the cc initial state in (b), for the db initial state in (c), and 
for the sb initial state in (d). The lower bounds on the new physics scale were taken from BHHS. 



improvement in the bounds could be found (or a discovery made) with relatively low integrated 
luminosity. Consider, for example, the uc initial state. There is currently no bound at all; 
in principle, A could be tens of GeV. The figure shows that a total integrated luminosity of 
an inverse picobarn would give a 5<r sensitivity for a A of 1 TeV. An integrated luminosity of 
an inverse femtobarn would give substantial improvements for all of the operators shown in 
Fig. |. 

4. Discussions and Conclusions 

In a previous article, motivated by discovery of large Vn — v T mixing in charged current 
interactions, bounds on the analogous mixing in neutral current interactions were explored. 
A general formalism for dimension-6 fermionic effective operators involving r — fi mixing with 
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typical Lorentz structure (jlTT)(q a Tq") was presented, and the low-energy constraints on 
the new physics scale associated with each operator were derived, mostly from experimental 
bounds on rare decays of r, hadrons or heavy quarks. Tensor operators were not considered, 
and some of the operators, such as cuflr, were completely unbounded. 

In this article, we consider \ir production at hadron colliders via these operators. Tables || 
and U list the new physics scales that are accessible at the Tevatron and the LHC, respectively. 
Due to much smaller backgrounds, both the LHC and Tevatron are more sensitive to electronic 
r decays than hadronic r decays. For hadronic r decays, the LHC receives an enhancement 
from the large invariant mass region and is more sensitive than the Tevatron. Since the 
backgrounds to electronic r decays at the Tevatron are much smaller than those at the LHC, 
the Tevatron is more sensitive than the LHC to electronic r decays. We found that at the 
Tevatron with 8 fb _1 , one can exceed current bounds for most operators, with most 2a 
sensitivities being in the 6 — 24 TeV range. We find that at the LHC with 1 fb _1 (100 fb _1 ) 
integrated luminosity, one can reach a 2a sensitivity for A ~ 3 — 10 TeV (A ~ 6 — 21 TeV), 
depending on the Lorentz structure of the operator. 
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A. New Physics Bounds 

The bounds from BHHS in units of TeV are presented in Table [|. The *s indicate there are 
no bounds on the new physics scale. Also, there are no bounds from BHHS for the tensor 
coupling. 



B. Partial Wave Unitarity Bounds 

Since the cross section from our higher-dimensional operators increases as s, it is necessary 
to determine the unitarity bound for qq — > [it. The partial wave expansion for a + b — > 1 + 2 
can be written as 

oo 

M(s,t) = 16vr ^2 (2J + l)aj(s)d^/(cos6>) 

J=M 

where 

1 f 1 

aj(s) = ^ J M(s, t)^^, (cos 9)d cos 9, 
[X = s a — Sb, [i' = si — S2 and J < max(|/i|, |//|). The condition for unitarity is |9i(aj)| < 1/2. 
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Coupling type 


1 


75 


7m 


7a*75 


uu 


2.6 


12 


12 


11 


dd 


2.6 


12 


12 


11 


ss 


1.5 


9.9 


14 


9.5 


ds 


2.3 


3.7 


13 


3.6 


db 


2.2 


9.3 


2.2 


8.2 


sb 


2.6 


2.8 


2.6 


2.5 


uc 


* 


* 


0.55 


0.55 


cc 


* 


* 


1.1 


1.1 


bb 


* 


* 


0.18 


* 



Table 5: Bounds on the new physics scales from BHHS in units of TeV for various operators and the 
scalar, pseudoscalar, vector, and axial-vector couplings. The *s indicate there were no bounds. 



It is straightforward to calculate the coefficients for the S,V,T operators. For example, 
for the scalar operator 

4-7T _ 

M = -pV Al (pi)uA 2 (p2)wA 3 (P3)wA4(P4) 



one can just plug in the explicit expressions: 



y /E - X\ P \ XX(P) 

^/E + X\p\xx(p) 



v\(p) = 



-y/E + \\p\ X -x(p)\ 

y/E- X\p\x-x(p) J 

where x+(z) = (o) iX-(^) = (?)• I n the massless limit, this simply gives ao = s/(4A 2 ) and 
so the unitarity bound gives s < 2A 2 . For the vector case, ao = and a\ = s/(6A 2 ) giving 
the unitarity bound s < 3A 2 . The tensor case gets contributions from both ao and a\, and 
the stronger bound then applies. 
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